Nuclear Structure Relevant to Neutrinoless Double Decay: 76 Ge and 76 Se. 
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The possibility of observing neutrinoless double /3 decay offers the opportunity of determining the 
neutrino mass if the nuclear matrix element were known. Theoretical calculations are uncertain and 
measurements of the occupations of valence orbits by nucleons active in the decay can be important. 
The occupation of valence neutron orbits in the ground states of 76 Ge and 76 Se were determined 
by precisely measuring cross sections for both neutron-adding and removing transfer reactions. Our 
results indicate that the Fermi surface is much more diffuse than in theoretical (QRPA) calculations. 
We find that the populations of at least three orbits change significantly between these two ground 
states while in the calculations the changes are confined primarily to one orbit. 

PACS numbers: 23.40.Hc, 25.40.Hs, 27.50.+e, 23.40-s 



An essential step in studying the nature of the neu- 
trino is the attempt to observe neutrinoless double beta 
decay [l[ and major efforts are being undertaken with this 
objective in mind. Observation of such a process would 
immediately show that neutrinos are their own antiparti- 
cles, and its rate may well give the first direct measure of 
the neutrino mass if the corresponding nuclear matrix el- 
ement can be reliably calculated. As an example, for one 
of the likely candidates ( 76 Ge), theoretical calculations 
have yielded answers that are spread over more than an 
order of magnitude. This prompted the statement by 
Bahcall et al. [3] "The uncertainty in the calculated nu- 
clear matrix elements for neutrinoless double beta decay 
will constitute the principal obstacle to answering some 
basic questions about neutrinos". There have been sug- 
gestions that relate the neutrinoless double-beta-decay 
matrix elements to those for ordinary single beta decay, 
or to the 'normal' two-neutrino modes which have been 
observed experimentally 0]. However, neutrinoless decay 
proceeds by the virtual excitation of states in the inter- 
mediate nucleus with a momentum transfer much larger 
than that for these other processes. It will thus involve 
all possible virtual intermediate states (up to about 100 
MeV of excitation), and so will include giant resonances. 
There is no other experimentally accessible process that 
could directly determine the matrix element. 

Although there is still considerable discussion regard- 
ing the best theoretical approach, what unquestionably 
matters is knowing the population of the valence orbits 
for the nucleons that switch from neutrons to protons. 
We have therefore undertaken a set of measurements to 
determine this quantity experimentally, and report here 
on the valence neutron populations and the differences in 
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these populations for 76 Ge and 76 Se. In a previous exper- 
iment we determined that the neutron pair correlations 
in these two nuclei are quantitatively very similar 4] . 

The Macfarlane- French [5( sum rules for nucleon trans- 
fer state that the summed spectroscopic strength for 
neutron- adding reactions with a given set of quantum 
numbers is equal to the vacancies in that target orbital, 
while the sum over states for neutron- removing reactions 
will determine the occupancy. Here we have measured 
the cross sections and extracted spectroscopic factors of 
significantly populated states, for both neutron-adding 
and neutron-removing reactions. The summed spectro- 
scopic factors for both reactions can be added and used 
to provide a normalization, allowing occupation numbers 
for orbit als to be extracted. 

The nucleon transfers reported here have been mea- 
sured previously [f| 0|, but not with the same experi- 
mental methods, and using different parameters in each 
DWBA analysis for extracting spectroscopic factors. The 
aim of the present measurement is to analyze all the re- 
sults in a consistent manner to permit the extraction of 
more accurate occupation numbers with a common ex- 
perimental approach. 

The active orbits for neutrons in these nuclei, with 
42 and 44 neutrons, are IP3/2, 0/ 5 / 2 , lpi/2i an d 0#9/2- 
We have made systematic measurements to obtain ac- 
curate cross sections for the neutron-adding (d,p) and 
(a, 3 He) reactions, as well as for neutron-removing (p,d) 
and ( 3 He,cv) reactions. The momentum matching in (d,p) 
reactions for transitions with £=3 and 4 is not optimal 
and thus the cross sections are rather weak. Therefore 
helium-induced reactions were used to obtain data with 
improved momentum matching and larger cross sections 
for the higher-^ transitions. This selectivity is illustrated 
in Figure 1. 

Deuteron, proton, alpha, and 3 He beams from the Yale 
tandem accelerator were used to bombard isotopically 
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FIG. 1: Energy spectra for the neutron-removal reactions for 
76 Se to 75 Se illustrating how different £ values are favored by 
different momentum transfers in the two reactions. The 1=1 
transitions appear strongly in the 11° (p,d) spectrum (points) 
while the £=3 and, in particular, ^=4 are most prominent in 
( 3 He,a) (line) where the resolution is worse because of the 
higher energy. The £ values are indicated in numbers above 
the peaks. 
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enriched Ge and Se targets of about 200-300 /xg/cm 2 
evaporated on thin, 50 ^g/cm 2 C foils. The momenta of 
the reaction products were determined and the particles 
identified with the Yale Enge spectrograph and gas-filled 
focal-plane detector backed by a scintillator. 

The product of target thickness and spectrometer solid 
angle was found by measuring elastic scattering in the 
Coulomb regime at 30° for each target used. The beam 
energies used for this calibration were 6-MeV protons 
and 10-MeV alphas. For the transfer reactions, the same 
spectrometer aperture and beam integrator settings were 
used to minimize potential systematic errors. The beam 
energies chosen were 15 MeV for the (d,p) reaction and 
23 MeV for the (p,d) to keep the energies in each channel 
comparable. Similarly, (a, 3 He) was studied at 40 MeV 
and ( 3 He,a) at 26 MeV. Measurements were also carried 
out on targets of 74 Ge and 78 Se to provide an additional 
check. The energy resolution obtained was ^40 keV for 
the deuteron and proton-induced reactions, and ~70 keV 
for the 3,4 He reactions. 



The (d,p) angular distributions have been studied pre- 
viously and £ values were assigned @, 0|- In the cur- 
rent work, the yields were therefore measured only at 
the angles that correspond to the peaks in the angular 
distributions for the t values of interest: 11°, 28° and 
37° for 1=1, 3 and 4 respectively. The helium- induced 
reactions are forward peaked and so the most practical 
forwardmost angles were chosen: 5° for (ev, 3 He) and 8° 
for its inverse. The previous £-value assignments d, 0] 
were confirmed, as may be seen in Figure 2. Our results 
also agree approximately with the previous relative spec- 
troscopic factors for states populated with a particular 
target. 



FIG. 2: Ratios of cross sections, a d>p (28°) /tJd, p (37°) 
vs. a d:P (ll )/a dtP (37°) on top and o a 3 Ht ,/a dtP {2^°) vs. 
°"d,p(ll )/°"ti,p(37 ) below, are shown for different £ values 
and reactions. The symbols, one for each state, indicate the 
^-value assignments from previous work: triangles (black) are 
1=1, circles (green online) are £=3, and stars (red online) are 
£=A. In addition, states not included in the analysis are £=2 
transitions indicated by x and £=0 by + signs. States with 
unknown I values are indicated by hollow circles (in blue on- 
line) . The size of the symbols is a rough measure of the cross 
sections. The dashed lines indicate the loci of the ratios for 
well-established £ values. The x surrounded by a (green on- 
line) circle, between the 1=2 and 3 islands in the lower box, 
is the 500-keV 5/2+ - 5/2" doublet in 77 Ge discussed in the 
text. 



We used the finite-range code PTOLEMY [| for the 
DWBA calculations. The normalization depends some- 
what on the choice of distorting parameters. The ex- 
tracted relative spectroscopic factors also vary slightly 
and this is a source of some of the uncertainty at the 
level of a few %. For the projectile bound-state wave 
function, the Reid potential was used for the deuteron, 
and a Woods-Saxon one for the a particle and for the 
various target bound states. 

Absolute spectroscopic factors are notoriously difficult 
to obtain. The values of spectroscopic factors for 'good' 
single-particle states in doubly-magic nuclei are usually 
around 0.5 - 0.6 because of short-range correlations. Such 
correlations are expected to be a uniform property of 
nuclei, not changing between nearby nuclei or configu- 
rations. Since the overall effect in depleting absolute 
strength is expected to be uniform, it can be compen- 
sated for by a renormalization of strength and examining 
the relative strengths of spectroscopic factors through the 
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sum rules. Since the sums of the strengths for neutron 
adding or removing are proportional to the vacancies or 
occupancies, together they should add up to the (2J+ 1) 
degeneracy of the orbits and can serve as such a normal- 
ization. A check is provided, in that the summed spectro- 
scopic factors for a given orbit should add up to the same 
value for each of the targets. The mean normalization 
factors with the potentials adopted for the 1=1, 3, and 4 
transitions are 0.53, 0.56, and 0.57 respectively with rms 
fluctuations among the targets of 2, 12, and 7 %, indicat- 
ing that the procedure is reasonable. This normalization 
constant for the two reactions is somewhat surprisingly 
close to the depletion that should be expected for 'abso- 
lute' spectroscopic factors. 

Several points are to be noted in the above sums. Since 
not all the spins of the states seen in £=1 transitions are 
known, we summed all 1=1 transitions, thus combining 
the j=l/2 and 3/2 states. For the neutron-removal reac- 
tions a small correction was made for the unobserved T> 
isobaric analog states, corresponding to proton removal. 
Also for neutron removal, in the (p,d) reaction, some 
previously determined £=1 transitions at high excitation 
energy were beyond the energy range measured here. A 
correction was made for these states by normalizing the 
previously determined spectroscopic factors to ones de- 
termined in the present work. There were no known 
missed states for the neutron-adding measurements or 
for the other £ values. 

Finally, for the f 5 / 2 states, no 5/2~ state was known in 
77 Ge, while all other nuclei in this region have such a state 
well below 1 MeV in excitation energy. In attempting to 
find such a state in the (a, 3 He) reaction, the intensity of 
the peak around 500 keV excitation was stronger than 
expected for a known £=2 transition to a 5/2 + state at 
504.8 keV, but the centroid of this peak seemed lower 
than the accepted value - around 492 keV. In fact, a 
tentative state is reported in the compilations @] at 491.9 
keV from unpublished work with the ( 13 C, 12 C) reaction, 
and we have assumed that this is the missing 5/2~ state. 
Its strength was included in the sums. 

The vacancies and occupancies from the summed nor- 
malized spectroscopic factors are shown in Table I. Listed 
in the Table are the numbers of holes and particles from 
neutron adding and removing, their sum, and the best 
average value of the occupancy, all computed with a con- 
stant normalization. The 1=1 strength is best deter- 
mined in the (d,p) and (p,d) reactions and the £=3 and 
4 transitions from the helium-induced reactions. As was 
noted, the sums of holes and particles for both £=1 and 
4 transfers are constant to better than 5% across the tar- 
gets studied. For £=3, the situation is somewhat worse, 
partly because these transitions are relatively weak in 
both reactions. As a result, components of strength could 
have been missed. Additionally, there is some ambigu- 
ity about the £=3 strength in the pickup reactions since 
some of the transitions could be to 7/2 ~ states at higher 
excitation energies. For cases where there is no evidence 
on the spins of higher £=3 hole states we, somewhat ar- 



TABLE I: Summed spectroscopic strengths. 
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bitrarily, excluded all £=3 transitions above 1 MeV ex- 
citation. As was noted, the summed strengths for £=3 
fluctuate more than the others. 

Our measurements provide two determinations of the 
valence-orbit occupancies in the 76 Ge and 76 Se ground 
states, one from the neutron-adding data, and one from 
the neutron removal. We average these, weighting the 
former by a factor of two. There are two reasons for this. 
First, the major neutron shell between N=28 and 50 is 
more than half filled, with about twice as many particles 
as holes. A given fractional error therefore leads to a 
bigger uncertainty in the particle number compared to 
the number of holes, at least for 1=1 and 3. Second, 
the sensitivity of the calculated DWBA cross sections to 
distorting parameters is larger for the ( 3 He,a) reaction 
(«15-25%) than for its inverse (w0-10%). 

The uncertainties in the final mean occupancy values 
are difficult to estimate. Statistical errors in the summed 
strength are less than 1% and relative systematic errors 
between targets are believed to be less than 3%. The 
biggest uncertainties stem from possible missed states, 
especially for the £=3 transitions, and from uncertainties 
in the DWBA calculations. We estimate that the occu- 
pancy is determined to about 0.2 nucleons for the lp, 0.3 
for the Ot/g/2 orbits, and slightly worse, 0.4, for the 0/ 5 / 2 
orbit. These estimates of uncertainties are rather crude. 
However, there are several checks that give us some con- 
fidence: 

• The normalization factors obtained for each target 
separately are similar. 

• The mean normalizations for each £ value and re- 
action type are also similar. 

• The summed removing and adding strengths for 
74 - 76 Ge and 76 ' 78 Se, 22.5, 20.7, 22.4, and 23.1 re- 
spectively, are consistent with the expected value 
of 22.0. 

• As an independent result, the neutron vacancies 
obtained for the four nuclei (from the adopted oc- 
cupancies in the Table) are 7.3, 6.1, 7.9 and 5.7, in 
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FIG. 3: The deduced neutron vacancies for 76 Ge and 76 Se 
are shown in the three active valence orbits and compared 
to those from the QRPA calculations of Reference 10. The 
naive shell closure should give 6 and 8 vacancies for these two 
nuclei. The lower part of the figure shows the differences in 
these occupations (expected to be 2.0), again compared to the 
QRPA calculation. 



good agreement with the expected values of 8, 6, 8, 
and 6. 

Beyond the valence lp, Of 5/2 and O59/2 orbitals, neu- 
tron removal from 76 Se suggests that approximately 0.2 
neutrons are in the 1^5/2 orbit. The weak 5/2 + state in 
75 Ge is not resolved in our work, but using the results in 



Q, we obtain a roughly similar value. 



The values of vacancies are shown in Figure 3 along 
with the QRPA results • There is little question that 
the vacancies in the lp and, especially, in the 0/ 5 / 2 or- 
bits are significantly larger in the data than in the cal- 
culations. For the neutrinoless double-beta-decay exper- 
iments it is the changes in occupancy that are important, 
and so in the lower part of Figure 3 we show the differ- 
ences between 76 Ge and 76 Se: 0.46±0.20 in lp, 0.73±0.40 
in 0/ 5/2 and 0.68±0.30 in 03 9/2 . 

While the QRPA results predict changes between the 
two nuclei to be mostly in the Ogg/2 orbit, the experiment 
shows quite clearly that the changes in the lp and 0/ 5 / 2 
orbits are much larger than predicted. The qualitative 
feature that, in disagreement with QRPA, there are still 
large vacancies in lp and O/5/2 is quite robust. It follows 
from the relatively large cross sections in the neutron- 
adding reactions and it cannot depend on the details of 
the analysis or the assumptions. 

What the consequences may be, of this disagreement 
in neutron occupancy between QRPA and experiment, 
on the matrix element for neutrinoless double beta decay 
are not clear at present and will need to be investigated in 
more detail. Proton occupancies are similarly important 
and experiments to determine them are planned. 
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sending us the results of the QRPA calculations. We 
wish to acknowledge the help of the operating staff of 
the Yale tandem and of John Greene at Argonne for 
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supported by the U.S. Department of Energy, Office of 
Nuclear Physics, under contracts DE-FG02-91ER-40609 
and DE-AC02-06CH11357, the UK Science and Technol- 
ogy Facilities Council, and IN2P3/CNRS-France. 
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